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Abstract — The  estimation  of  temperature  inside  battery  cells 
requires  accurate  information  about  the  cooling  conditions  even 
when  the  temperature  of  the  battery  surface  is  measured.  This 
paper  presents  a  novel  approach  of  estimating  temperature  dis¬ 
tribution  inside  cylindrical  batteries  under  unknown  convective 
cooling  conditions.  A  computationally  efficient  thermal  model 
is  first  developed  using  a  polynomial  approximation  of  the 
temperature  profile  inside  the  battery  cell.  The  Dual  Extended 
Kalman  Filter  (DEKF)  is  then  applied  for  the  identification  of 
the  convection  coefficient  and  the  estimation  of  temperature 
inside  the  battery.  In  the  proposed  modeling  approach,  the 
thermal  properties  are  represented  by  volume  averaged  lumped 
values  with  uniformly  distributed  heat  generation.  The  model 
is  parameterized  and  validated  using  experimental  data  from  a 
2.3  Ah  26650  Lithium-Iron-Phosphate  (LFP)  battery  cell  with  a 
forced-air  convective  cooling  during  hybrid  electric  vehicle  (HEV) 
drive  cycles.  Experimental  results  show  that  the  proposed  DEKF- 
based  estimation  method  can  provide  an  accurate  prediction  of 
core  temperature  under  unknown  cooling  condition  by  measuring 
the  cell  current,  voltage,  and  surface  and  ambient  temperature. 
The  accuracy  is  such  that  the  scheme  cam  also  be  used  for  fault 
detection  of  a  cooling  system  malfunction. 

Index  Terms — Lithium  ion  batteries,  Thermal  model,  Reduced 
order  model,  Dual  Extended  Kalman  Filter,  Parameter  identifi¬ 
cation 

I.  INTRODUCTION 

OVER  the  past  years,  energy  storage  systems  utilizing 
lithium  ion  (Li-ion)  batteries  have  become  one  of  the 
most  critical  components  for  realizing  efficient  and  clean 
transportation  systems  through  electrification  of  vehicles,  e.g., 
hybrid  electric  vehicles  (HE Vs),  plug-in  hybrid  electric  ve¬ 
hicles  (PHEVs),  and  electric  vehicles(EVs).  Li-ion  batteries 
have  several  advantages  -  no  memory  effect,  wide  range  of 
operating  temperature,  and  high  energy  and  power  density 
[1],  [2].  However,  the  Li-ion  battery  performance,  cycle  life 
and  capacity  are  adversely  affected  by  sustained  operation  at 
extreme  (above  45°C  and  below  freezing)  temperatures  [3]- 
[6],  a  recurring  problem  in  automotive  applications  where 
batteries  are  exposed  to  temperature  extremes  with  frequent 
high  current  discharge/charge  rate  that  cause  internal  heating. 
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Thus,  being  able  to  estimate/predict  the  temperature  distri¬ 
bution  across  cells  and  packs  is  vital  for  formulating  power 
management  strategies  that  are  mindful  of  the  performance 
limitations  of  these  versatile  power/energy  sources.  In  general, 
the  performance  of  the  cooling  system  can  be  degraded  due 
to  various  reasons  such  as  dust  on  fan  blades,  partial  blockage 
in  pipes,  motor/pump  ageing,  and  even  a  motor/pump  failure. 
When  such  a  degradation  or  failure  occurs,  it  is  not  possible 
to  reject  the  heat  generated  from  the  battery  cell.  In  this 
condition,  the  lifespan  of  the  battery  exposed  to  the  extreme 
temperatures  will  be  considerably  shortened.  Therefore,  it  is 
important  to  identify  the  convective  heat  coefficient  not  only 
to  accurately  estimate  the  temperature  distribution  inside  the 
battery  but  also  for  fault  detection  to  ensure  safe  and  reliable 
operation  of  the  vehicle  system. 

This  paper  considers  a  novel  method  for  estimating  temper¬ 
ature  distribution  inside  cylindrical  batteries  with  simultaneous 
estimation  of  the  convective  cooling  condition.  To  achieve 
this  goal,  a  computationally  efficient  thermal  model  for  a 
cylindrical  cell  is  utilized  to  estimate  simultaneously  the 
convection  coefficient  and  radial  temperature  distribution  of 
the  cell.  Unlike  existing  reduced  order  modeling  approaches 
in  [7] — [10],  and  [11],  a  polynomial  approximation  to  the 
solution  of  the  heat  transfer  problem  is  used;  this  approach 
facilitates  a  systematic  estimation  of  core,  surface,  volume- 
averaged  temperatures,  and  volume-averaged  temperature  gra¬ 
dients.  Dual  Extended  Kalman  Filter  (DEKF)  is  applied  for  the 
identification  of  the  convection  coefficient  and  the  temperature 
distribution  inside  a  cylindrical  battery  cell.  The  proposed 
estimation  method  provides  the  capability  of  detecting  the 
malfunction  of  cooling  system  by  monitoring  the  difference 
between  the  identified  and  off-line  predetermined  convection 
coefficient.  This  benefit  indicates  that  a  significant  rise  in 
temperature  can  be  prevented  by  augmenting  the  proposed 
method  with  other  existing  battery  management  strategies  for 
the  safe  and  robust  operation. 

This  paper  is  organized  as  follows:  Section  II  presents  the 
convective  heat  transfer  problem  for  a  cylindrical  battery  cell 
and  the  reduced  order  thermal  model.  Model  reduction  is 
performed  using  a  polynomial  approximation  of  the  partial 
differential  equation  (PDE)  system.  The  thermal  properties 
of  the  battery  are  experimentally  identified  and  sensitivity 
of  parameters  is  analyzed  in  Section  III.  In  Section  IV, 
the  temperature  estimator  applying  a  Dual  Extended  Kalman 
filter  by  using  the  proposed  model  for  estimating  the  core 
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Fig.  1.  Schematic  for  a  A 123  26650  cylindrical  battery  cell 


temperature  and  identifying  the  convection  coefficient.  Section 
V  presents  and  discusses  experimental  results  and  conclusions 
are  drawn  in  Section  VI. 


II.  HEAT  TRANSFER  PROBLEM  IN  CYLINDRICAL 
BATTERIES 


This  paper  considers  the  radially  distributed  (1-D)  thermal 
behavior  of  a  cylindrical  battery  cell  with  convective  heat 
transfer  boundary  condition  as  illustrated  in  Fig.  1  [7],  [12], 
[13],  A  cylindrical  Li-ion  battery,  so-called  a  jelly-roll,  is  fab¬ 
ricated  by  rolling  a  stack  of  cathode/separtor/anode  layers.  The 
individual  layered  sheets  are  thin,  therefore,  it  is  reasonable 
to  assume  uniform  heat  generation  along  the  radial  direction 
[13],  [14].  Lumped  parameters  are  used  so  that  material 
properties  such  as  thermal  conductivity,  density,  and  specific 
heat  coefficient  are  assumed  to  be  constant  in  a  homogeneous 
and  isotropic  body.  The  thermal  conductivity  is  one  or  two 
orders  of  magnitude  higher  in  the  axial  direction  than  in 
the  radial  direction.  Therefore,  the  temperature  distribution 
in  the  axial  direction  will  be  more  uniform  [15],  [16].  The 
governing  equation  of  the  1-D  temperature  distribution  T(r,  f) 
and  boundary  conditions  are  given  by 


dT(r,t )  d2T(r,t )  kthdT(r,t )  Q(t) 

Pcp — 777 — =  kth — - 1 - 77 - 1 — ^ — ,  (1) 


dt 

dT(r,  t) 
B.C.’s  — , '  7 
or 

dT(r,t ) 


dr2 


dr 


Vb 


=  0, 


r— 0 


(2) 


dr 


=  -  —  (T(R,t)-T00),  (3) 

r — R  kth 


where  t,  p.  cp  and  kth  represent  time,  volume-averaged  den¬ 
sity,  specific  heat  coefficient,  and  thermal  conductivity  of  the 
cell  respectively.  The  radius  of  the  battery  cell  is  R,  Q  is 
the  heat  generation  inside  the  cell,  and  VJ,  is  the  volume  of 
battery  cell.  Ambient  temperature  for  convection  is  denoted  by 
Too.  The  boundary  condition  in  (2)  represent  the  symmetric 
structure  of  the  battery  about  the  core.  The  other  boundary 
condition  shown  in  (3)  represents  the  convective  heat  transfer 
at  the  surface  of  the  battery. 


A.  Model  reduction 

With  evenly  distributed  heat  generation,  the  temperature 
distribution  along  r-direction  of  the  battery  cell  is  assumed 


to  satisfy  the  following  polynomial  approximation  proposed 
in  [17] 

T(r,  f)  =  a(f)  +  6(f)  +dW(^)  >  (4) 

where  a(t),  6(f),  and  d(t)  are  time-varying  constants.  To 
satisfy  the  symmetric  boundary  condition  at  the  core  of 
the  battery  cell,  (4)  contains  only  even  powers  of  r.  Thus, 
the  temperatures  at  core  and  surface  of  the  battery  can  be 
expressed  as 

T(0,f)  =Tc=a(f),  (5) 

T(R,t)  =TS  =a(t)  +  b(t)  +  d(t),  (6) 


where  subscripts  c  and  s  denote  core  and  surface  respectively. 

The  volume-averaged  temperature  T  and  temperature  gra¬ 
dient  7  are  introduced  as  follows: 

2  rR 

T  =  R2lrTdr ’  (?) 


7  = 


2  fR  ( dT\ 

&  Jo  r\fr) 


dr. 


(8) 


These  volume-averaged  values  are  used  as  the  states  unlike 
existing  approaches  in  [8],  [9],  and  [11]. 

By  substituting  (4)  in  (7)  and  (8),  T  and  7  can  be  expressed 
in  terms  of  constants  as 


T  =  a(t)+bM  +®, 

46(f)  ^  8 d(t) 


3  R 


5  R 


(9) 

(10) 


By  rearranging  (6),  (9),  and  (10),  time-varying  constants 
a(f),  6(f),  and  d(t)  can  be  written  by 

a(f)=4Ts-3T - —7,  (11) 

15  P 

6(f)  =  -  18TS  +  18T  +  —7,  (12) 

d(t)  =  15TS  -  15T  -  ^7.  (13) 

O 

By  substituting  (11),  (12),  and  (13)  in  (4),  the  temperature 
distribution  can  be  expressed  as  a  function  of  Ts,  T.  and  7 


1 

T(r,  t)  =  4TS  —  3T - 


-18Ts  +  18T+^7 

,r,r  ,rr-r  45 RJ  ' 
15  ls  —  15i - —7 


(s)! 


(- 

\R 


(14) 


The  PDE  (1)  can  be  converted  into  ODEs  by  substituting 
(14)  in  volume-averaged  governing  equation  and  its  partial 
derivative  with  respect  to  r  as  follows: 


dT 

dt 


48a - 
R2 


48a 

R2 


Ts 


dj  320a  - 
dt  R3 


15  a. 
"  R  ' 
320a 
R3 


a 


<2  =  0, 


kthVb 
120  a_ 


(15) 

(16) 


where  a  is  thermal  diffusivity  and  is  defined  as  follows: 


a  —  kth/  pCp. 


(17) 
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Fig.  2.  Comparison  of  frequency  response  functions  between  analytical 
solution  and  polynomial  approximation 


Using  (3),  the  surface  temperature  Ts  can  be  rewritten  as 

T  =  24  kth  j,  15  kthR  _  Rh 

s  24 kth  +  Rh  48 kth  +  ZRh1  2Akth  +  Rh  °°' 

08) 

Finally,  a  two-state  thermal  model  can  be  given  by  the 
following  form: 


x  =  Ax  +  Bu, 

y=Cx  +  Du,  (19) 


where  x  =  [T  ^]T ,  u  =  [Q  T^Y  and  y  =  [Tc  Ts]t  are 
states,  inputs  and  outputs  respectively.  System  matrices  A,  B, 
C,  and  1)  are  defined  as  follows: 


A  = 

B  = 

C  = 

D  = 


—48  ah  —15  ah 

R(24kth~\~Rh)  24kth~\~Rh 

-320  ah  -120a(4kth+Rh)  ) 

R2(24kth+Rh)  R2{24kth  +  Rh) 

cx.  _ 48  ah _ 

kth  Fj)  R.(24kth~\-Rh) 

n  _ 320a: /t _  ? 

U  R2(24kth+Rh) 

24kth-3Rh  120Rkth  +  15R2h 
24kth4~Rh  8(24  kth~\~Rh) 

24  kth  15  Rkth  > 

24  kth+Rh,  48kth+2Rh 

rj  4  Rh 

24k±h+Rh 
n  Rh 

u  24kth+Rh  . 


(20) 


This  state-space  representation  is  used  for  the  development  of 
control  design. 


TABLE  I 

Parameters  of  the  battery  [7] 


Parameter 

Symbol 

Value 

Unit 

Density 

P 

1824 

kg/m3 

Specific  heat  coeff. 

Cp 

825 

J/kgK 

Thermal  conductivity 

kth 

0.488 

W/mK 

Convection  coeff. 

h 

5 

W/m2-K 

Radius 

R 

12.93e-3 

m 

Height 

L 

65.15e-3 

m 

Volume 

vb 

3.4219e-5 

m3 

rameters  used  to  generate  the  plots  in  Fig.  2  are  summarized  in 
Table  I.  The  heat  transfer  coefficient  of  /i=5W/m2K  is  chosen 
since  this  value  is  typical  of  natural  convection  condition  [18]. 

Figure  2  shows  that  the  effects  of  heat  generation  on  core 
and  surface  temperature,  denoted  by  H\  -\  (s)  and  ll^t  (a)  re¬ 
spectively,  can  be  accurately  predicted  over  the  whole  range  of 
frequency.  On  the  other  hand,  the  responses  of  core  and  surface 
temperature  excited  by  the  ambient  temperature,  //|  2  (a)  and 
H22(s),  are  nearly  identical  to  the  analytical  solution  for 
frequencies  below  10-2  Hz.  In  general,  the  temperature  of 
cooling  media  does  not  change  rapidly;  thus,  the  prediction  of 
temperature  distribution  using  the  proposed  approach  can  be 
considered  sufficiently  accurate. 

C.  Heat  generation  calculation 

Since  heat  generation  rate  Q  is  the  input  to  the  battery 
thermal  system,  the  input  needs  to  be  accurately  calculated 
from  measurement  data,  such  as  current  and  voltage  during 
operation.  In  [19],  Bernardi  et  al.  proposed  the  simplified  form 
of  heat  generation  rate  with  assumptions  that  heat  generation 
due  to  enthalpy-of-mixing,  phase-change,  and  heat  capacity 
are  assumed  to  be  negligible  expressed  by 

Q  =  i(u  -  V)  -  i  (r?pj  ,  (22) 

where  i,  U,  and  V  represent  the  current,  the  open-circuit 
voltage  (OCV),  and  the  terminal  voltage  respectively.  As 
shown  in  Fig  3,  the  OCV  is  a  function  of  the  battery  state- 
of-charge(SOC).  This  function  is  experimentally  obtained  by 
averaging  the  measured  terminal  voltages  during  charging 
and  discharging  a  battery  with  C/20  current  rate  under  a 
Constant  Current  Constant  Voltage  (CCCV)  charging  protocol. 


B.  Frequency  domain  analysis 

The  transfer  function  of  the  thermal  system  II (s')  is  calcu¬ 
lated  by 

H(s)  =D  +  C(sI-  A)~1B,  (21) 

where  s  is  Laplace  variable. 

The  frequency  response  of  transfer  function  of  the  proposed 
model  is  compared  to  that  of  the  analytical  solution  in  [7].  Pa¬ 


Fig.  3.  Open-circuit  voltage  approximately  obtained  by  averaging  terminal 
voltages  during  charging  and  discharging  a  battery  with  C/20  current  rate 
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Fig.  4.  Data  set  used  for  parameter  ID:  current  (top),  voltage  (middle),  and 
heat  generation  rate  (bottom)  during  Urban-Assault  Cycle 


The  OCV  is  then  calculated  at  the  estimated  SOC  value  by 
integrating  measured  current  with  respect  to  time  as 


dSOC  I 

dt  3600Qh 


(23) 


where  Qi,  is  a  battery  capacity  in  Ah.  The  sign  convention  is 
such  that  positive  current  denotes  battery  discharging. 

The  last  term  in  (22)  is  the  heat  generation  from  entropy 
change.  In  this  paper,  heat  generation  due  to  entropy  change 
is  neglected  for  simplicity.  This  simplification  is  warranted 
since  the  typical  SOC  range  of  HEV  operation  is  between 
40%  and  60%  in  which  ^  of  the  battery  cell  is  insignificant 
as  shown  in  [9]  for  this  chemistry.  In  addition,  the  reversible 
entropic  heat  generation  would  have  zero  mean  value  when 
the  battery  is  operating  in  charge-sustaining  mode,  typical  of 
HEV  operation. 


III.  PARAMETER  IDENTIFICATION 

In  this  section,  the  value  of  the  lumped  parameters  in  (19) 
for  a  2.3  Ah  26650  LFP  battery  cell  by  A123  are  identified 
through  experimentation.  Figure  4  shows  current,  voltage  and 
calculated  heat  generation  rate  profiles  over  power  demanding 


TABLE  II 

Identified  thermal  properties 


Parameter 

Symbol 

Value 

Reference 

Density 

P 

2047* 

2118  [20] 

Sp.  heat  coeff. 

Cp 

1148.1 

1004.9-1102.6  [9],  [21] 

Thermal  cond. 

kth 

0.698 

0.69  [20] 

Conv.  coeff. 

h 

60.00 

65.99  [21] 

*  Measured 


Time(sec) 

Fig.  5.  Comparison  between  measured  and  simulated  temperatures  (top) 
and  errors  (bottom) 


cycle  used  for  military  ground  vehicle  [22]  Urban-Assault 
Cycle  (UAC)  that  are  used  for  the  parametrization.  The  model 
is  then  validated  using  a  different  duty  cycle.  The  numerical 
analysis  on  parameter  sensitivity  is  performed  to  investigate 
the  use  of  constant  parameters  for  thermal  conductivity  and 
heat  capacity  and  the  importance  of  identifying  the  convection 
coefficient  on-line. 

A.  Identifying  thermal  properties 

Parameter  identification  is  important  for  accurately  pre¬ 
dicting  the  temperature  distribution  inside  a  battery  cell  as 
the  parameters  kth ,  cp,  p,  and  h  determine  the  dynamics 
of  thermal  model.  As  the  density  can  be  assumed  to  be 
a  measurable  constant,  only  three  parameters  such  as  kth , 
cp,  and  h  are  considered  for  the  parameter  identification. 
Following  an  experimental  set-up  in  [21],  we  measured  a  cell 
current,  voltage,  surface  and  core  temperature  of  the  battery 
cell,  and  ambient  temperature  for  parameter  identification.  As 
discussed  in  [21],  some  parameters  are  assumed  to  be  known 
since  all  parameters  are  not  identifiable  without  using  the  core 
temperature. 

Let  the  error  between  the  measured  temperatures  and  model 
outputs  at  each  time  step  k  in  vector  form  be 

e(M)  =  [Tc,e(k,9)  TSte(k,0)}T  -  [Tc,m(fc)  Ts,m{k))T , 

(24) 

where  6  =  [kth  cp  h]T,  TC:B  and  Ts  e  represent  the  model 
parameters,  core  and  surface  temperatures  respectively. 

Parameters  are  identified  by  minimizing  the  Euclidean  norm 
of  the  difference  between  the  measured  and  simulated  temper¬ 
atures  as  given  by 

Nf 

9*  =  arg  min  E  ||e(M)||2,  (25) 

k= 1 

where  Nf  is  the  number  of  measurement  points.  The  mini¬ 
mization  problem  is  solved  by  using  the  finincon  function  in 
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Fig.  6.  (a)  Expected  temperature  distribution  along  the  normalized  radius  ( r/R )  using  polynomial  approximation;  (b)  Cell  temperature  (top)  and  temperature 

gradient  (bottom) 


MATLAB;  the  parameters  in  Table  I  are  used  as  initial  guess 
for  the  identification. 

Table  II  presents  the  identified  thermal  properties  for  the 
A123  26650  battery;  these  parameters  are  close  to  the  values 
presented  in  the  literature.  The  identified  specific  heat  coeffi¬ 
cient  cp  is  4%  larger  than  the  maximum  value  determined  in 
[9]  where  cp  was  determined  by  measuring  transient  responses 
of  the  battery  under  different  pulses.  It  was  discussed  that 
the  deviation  in  identified  value  of  cp  might  be  caused  by 
measurement  uncertainty  in  temperature  and  the  temperature 
dependency  of  the  heat  capacity. 

Despite  using  similar  experimental  data  and  setup,  the  iden¬ 
tified  convection  coefficient  is  10%  smaller  than  the  coefficient 
calculated  by  using  thermal  resistance  and  battery  surface 
area  in  [21].  This  difference  between  our  identified  value 
and  the  one  in  [21]  may  be  due  to  the  two  different  model 
structures.  Lin  et  al.  in  [21]  considered  two  different  materials, 
namely  one  for  the  core  and  one  for  the  surface,  whereas 
we  assume  the  battery  is  a  homogeneous  and  isotropic  body. 
In  order  to  accurately  determine  the  convection  coefficient, 
the  temperature  measurements  of  a  pure  metal  during  thermal 
relaxation  can  be  used.  For  instance,  the  specific  heat  capacity 
of  copper  at  25°C  is  known  as  385  J/kgK.  For  more  detailed 
description  about  the  experiment,  the  interested  reader  is 
referred  to  [23], 

Figure  5  shows  the  measured  and  simulated  temperatures 
at  the  core  and  surface  of  the  battery.  The  error  between 
the  measurements  and  simulated  temperature  is  less  than  the 
sensor  accuracy  of  0.5°C.  Thermocouples  used  for  temperature 
measurements  are  T-type  whose  accuracy  is  the  maximum 
of  0.5°C  or  0.4%  according  to  technical  information  from 
the  manufacture,  OMEGA.  The  parameterized  thermal  model 
accurately  predicts  the  temperature  inside  the  battery,  which  is 
difficult  to  measure  in  practical  applications.  Using  (15),  the 
temperature  distribution  inside  the  battery  can  be  predicted  as 
presented  in  Fig  6(a).  Figure  6(b)  shows  the  volume-averaged 
temperature  and  its  gradient  of  the  battery  respectively.  As 


evidenced  in  Fig.  6(b),  there  is  no  significant  difference  be¬ 
tween  the  volume-averaged  temperature  and  the  linear  average 
of  core  and  surface  temperatures,  i.e.  ( Ts  +  Tc)/2.  It  should 
be  noted  that  existing  approaches  in  [8],  [9],  [11]  have  the 
capability  of  predicting  the  core  temperature  and  have  shown 
the  efficacy  of  their  proposed  methods  on  the  prediction  of 
temperature  inside  the  cell  under  consideration  in  this  work. 
However,  the  phenomena  may  differ  in  the  case  of  a  cell 
with  larger  radius  [24],  The  volume-averaged  temperature 
gradient  is  different  from  the  linear  temperature  gradient,  i.e. 
(Ts  —  Tc)/R\  in  particular,  the  volume-averaged  temperature 
gradient  is  1.36  times  greater  than  linear  temperature  gradient 
under  the  UAC  test.  Since  non-uniform  temperature  distri¬ 
bution  can  lead  to  accelerated  capacity  losses  of  inner  core 
[24],  the  volume-averaged  temperature  gradient  is  an  important 
metric  to  describe  severity  of  temperature  inhomogeneity 
inside  the  battery. 

B.  Model  validation 

In  order  to  validate  the  performance  of  the  proposed  model 
with  identified  parameters,  the  battery  was  tested  under  a 
different  HEV  drive  cycle,  the  Escort  Convoy  Cycle  (ECC). 
The  current  and  voltage  profiles  for  this  cycle  are  illustrated 
in  Fig  7(a).  Figure  7(b)  shows  that  there  are  slight  differences 
between  the  measured  and  simulated  temperatures;  in  partic¬ 
ular,  the  root-mean-square  errors  (RMSE)  of  core  and  surface 
temperatures  are  0.4  and  0.3°C  respectively.  These  differences 
may  be  explained  with  the  assumption  of  radially  uniform 
heat  generation  and  high  conductivity  in  the  axial  direction. 
Additionally,  the  hysteresis  effect  of  the  LFP  battery  is  not 
properly  considered  in  heat  generation  formulation  (22),  which 
might  introduce  error  in  the  calculation  of  heat  generation 
rate.  Nevertheless,  since  the  comparison  of  temperatures  shows 
good  agreement  and  reasonably  small  RMSEs,  it  can  be 
concluded  that  the  proposed  model  with  identified  thermal 
properties  is  sufficiently  accurate  for  HEV  drive  cycles. 
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Fig.  7.  (a)  Current  (top)  and  voltage  (bottom)  profiles  during  Escort  Convoy  cycle;  (b)  Validation  data  set:  comparison  between  measured  and  simulated 

temperatures  (top)  and  errors  (bottom) 


C.  Parameter  sensitivity  analysis 

In  order  to  investigate  the  impact  of  parameter  variations 
on  the  performance  of  temperature  prediction,  each  parameter 
is  varied  from  the  identified  value  while  holding  the  other 
parameters  constant.  Figure  8  shows  that  parameters  such 
as  thermal  conductivity  kth  and  specific  heat  capacity  cp 
have  more  influence  on  the  prediction  of  core  temperature 
than  surface  temperature.  This  result  corresponds  to  the  fact 
that  the  heat  inside  the  battery  cell  is  transferred  through 
the  conduction.  On  the  other  hand,  the  prediction  of  surface 
temperature  is  most  sensitive  to  the  variation  of  convection 
coefficient,  which  can  be  explained  given  the  fact  that  the 
convection  coefficient  is  directly  related  to  the  convective 
boundary  condition  (3).  The  convection  coefficient  has  the 
most  significant  influence  on  the  overall  prediction  of  the  core 
and  surface  temperature. 

The  specific  heat  coefficient  and  thermal  conductivity  are 


Variation  (%) 

Fig.  8.  The  effect  of  parameter  variation  to  the  prediction  of  core  and  surface 
temperatures 


weakly  dependent  on  temperature  [9],  [25],  [26],  so  the 
assumption  of  constant  parameters  can  be  justified.  On  the 
other  hand,  the  convection  coefficient  is  highly  dependent 
on  fan  speed  or  fluid  velocity  as  expressed  by  empirical 
correlations  provided  by  Zukauskas  [27].  Consequently,  the 
accurate  identification  of  convection  coefficient  is  important 
for  better  prediction  of  temperature  inside  the  battery.  This 
importance  justifies  the  on-line  identification  of  the  convection 
coefficient  for  better  estimation  of  temperature  as  detailed  in 
Section  IV. 


IV.  ESTIMATION  OF  TEMPERATURE  AND 
CONVECTION  COEFFICIENT 

As  discussed  in  section  III-C,  the  estimation  of  temperature 
inside  the  battery  cell  requires  accurate  knowledge  of  the 
convection  coefficient  which  depends  on  cooling  condition. 
In  order  to  identify  the  convection  coefficient  on-line,  the 
Dual  Extended  Kalman  filter  (DEKF)  [28]  is  applied  for  better 
estimation  of  temperature  distribution  inside  the  battery  cell. 
The  other  thermal  parameters,  such  as  thermal  conductivity 
and  specific  heat  coefficient,  have  less  impact  on  temperature 
and  do  not  change  significantly  over  time.  Therefore,  the 
constant  values  identified  in  section  III  can  be  used. 

Assuming  the  input  u(t)  is  constant  over  each  sampling 
interval,  a  parameter  varying  (PV)  discrete-time  model  at  time 
step  k  can  be  obtained  as 

‘t'fc+l  ^di.9k)Xk  T  Bd{@k)'U'ki 

Uk  =  C{9k)xk  +  D(9k)uk,  (26) 

where  x  =  [T  j]7 ,  y  =  [Tc  Ts j1 ,  9  =  h,  and  u  =  [Q  Too]71. 
System  matrices  Ad  ~  I+AAT  and  Bd  =  BAT  are  obtained 
from  matrices  in  (20)  where  the  sampling  period  is  AT,  and 
I  is  the  identity  matrix. 

Let  the  PV  thermal  system  in  discrete-time  domain  be 
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Fig.  9.  Fan  schedule  for  forced- air  convective  cooling 


expressed  in  a  general  form  by: 

^k-\- 1  f  @k)  "b  ^ki 

Vk  =  g{xk,uk,dk)  +  vk , 

0k+i  =  0k  +  rk,  (27) 


where  wk,  vk,  and  rk,  are  independent,  zero-mean,  Gaussian 
noise  processes  of  covariance  matrices  Ew,  E„,  and  Er, 
respectively.  The  design  of  the  DEKF  estimator  is  given  as 
following  update  processes. 

Time  update  for  the  parameter  filter: 

ez=ei_i1  (28) 

Sk  =  S+_ !  +  Er.  (29) 

Time  update  for  the  state  filter: 

=  f(xi-i,y-k-udk)>  (30) 

P-  =Ak-1P+_1Al_1  +  Zw.  (31) 


Measurement  update  for  the  state  filter: 


Kk=P-Cf 


=*fc 


K, 


CIPkCf  +  E 
Vk  - a(xk,uk,ek) 


Pk  =  [I  ~  KkCk  \  Pk  ■ 


(32) 

(33) 

(34) 


Measurement  update  for  the  parameter  filter: 


Lk  =  Skcf 


ckPkckT  + 


I  “I 


Lk 


Vk-  g(xk,uk,ok ) 


[l-LkCek\Sk, 


(35) 

(36) 

(37) 


where  the  matrices  are  calculated  by 

_df(xk-i,uk-1,ek) 
k~1  -  dxk-1  lxfe-i=* 
Cx  =  dg{xk,uk,6k) 
k  dxk  lxf==£fc  ’ 

r,e  _  dg(xk,uk,d)  , 

“  d6  le=^_  ' 


+ 

fc-1 


(38) 

(39) 

(40) 


The  details  of  the  matrices  are  provided  in  Appendix  A. 
Superscripts  —  and  +  denote  the  a  priori  and  a  posteriori 
values  respectively. 

The  identified  states  x  and  parameter  8,  computed  from  the 
above  DEKF  algorithm,  are  used  to  estimate  the  core  temper¬ 
ature  in  the  battery  cell  from  (26).  The  identified  parameter 
can  be  also  used  for  monitoring  the  malfunction  or  degradation 


Fig.  10.  Comparison  of  performance  between  KF  estimator,  DEKF  estimator 
during  stage  I:  convection  coefficient  (top)  and  temperature  (bottom) 


of  cooling  system.  Under  the  assumption  that  the  relationship 
between  the  convection  coefficient  and  fan  speed  or  PWM 
signal  is  known,  the  malfunction  of  the  cooling  system  can 
be  detected  by  comparing  the  identified  parameter  with  a 
predetermined  range  of  values  8*  based  on  the  identification 
process  performed  at  various  cooling  conditions.  When  the 
difference  between  the  identified  and  predetermined  values 
\6  —  8*  |  is  bounded  and  small,  it  can  be  considered  that  there 
is  no  fault  in  the  cooling  system.  On  the  other  hand,  where 
\8  —  8*\  3>  e  where  e  is  a  pretuned  threshold,  a  cooling  fault 
can  be  detected.  Furthermore,  \{0  —  0*)/6*\  can  be  interpreted 
as  the  severity  of  degradation  of  cooling  system. 

V.  EXPERIMENTAL  RESULTS 

In  this  section,  the  performance  of  the  proposed  temperature 
estimator  using  the  DEKF  is  compared  with  that  of  the 
baseline  Kalman  Filter  (KF)  estimator  without  parameter  iden¬ 
tification.  The  battery  is  tested  using  the  ECC  under  different 
cooling  conditions.  Three  different  forced  convective  cooling 
conditions  (stage  I,  stage  II,  and  stage  III)  are  achieved  by 
using  different  PWM  signals  driving  the  fan  as  shown  in  Fig. 
9  which  corresponds  to  an  increase,  followed  by  a  decrease 
in  the  coolant  flow  rate.  In  order  to  investigate  the  influence 
of  change  in  the  parameter  on  the  temperature  estimation,  the 
parameter  is  provided  to  each  estimator  as  following: 

•  In  stage  I,  the  off-line  predetermined  convection  coeffi¬ 
cient  is  provided  to  the  KF  and  is  used  for  the  DEKF  as 
an  initial  value:  8  =  8*  and  0(0)  =  8* 

•  In  stage  II,  the  off-line  predetermined  convection  coeffi¬ 
cient  is  provided  to  the  KF  only:  9  =  8* 

•  In  stage  III,  two  times  larger  convection  coefficient  com¬ 
pared  to  the  known  value  is  provided  to  the  KF:  8  =  29* 

where  9  and  8  denote  fixed  and  identified  parameters  for 
the  KF  and  the  DEKF  respectively.  Other  thermal  properties 
such  as  thermal  conductivity  and  specific  heat  coefficient  are 
assumed  constant  with  values  identified  in  section  III. 
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Fig.  11.  Comparison  of  performance  between  KF  estimator,  DEKF  estimator 
during  stage  II:  convection  coefficient  (top)  and  temperature  (bottom) 


5000  5500  6000  6500  7000 


Time(sec) 

Fig.  12.  Comparison  of  performance  between  KF  estimator,  DEKF  estimator 
during  stage  III:  convection  coefficient  (top)  and  temperature  (bottom) 


It  is  assumed  that  the  initial  temperature  distribution  inside 
the  battery  is  uniform  at  30°C  and  convection  coefficient  is 
56.2  W/mn2K,  i.e.  £(0)  =  [30  0]7  and  0(0)  =  56.2  respec¬ 
tively.  The  covariance  matrix  for  the  state  Ew  =  /?i2diag(l,  1) 
describes  the  process  noise  where  /3  >  0  is  a  parameter 
for  tuning  relaxed  to  model  inaccuracy.  The  noise  covariance 
E„  =  a2  is  determined  from  the  standard  deviation  of 
temperature  signal  a  =  0.05°C.  The  covariance  matrix  for 
the  parameter  Er  =  /322  influences  the  performance  of  noise 
filtering  and  the  rate  of  parameter  convergence.  Ultimately,  the 
initial  condition  of  the  error  covariance  matrix  and  the  tuning 
parameter  are  chosen  as  P( 0)  =  diag(l,l),  /3i  =  0.0005, 
S'(O)  =  1,  and  /32  =  0.01  through  repeated  simulations 
respectively. 

The  results  for  the  parameter  and  state  estimation  are 
shown  in  Fig.  10-12  and  summarized  in  Table  III.  Figure  10 
shows  that  the  closed  loop  estimators  can  accurately  predict 
temperature  inside  the  battery.  Even  though  the  identified 
value  is  used  as  an  initial  guess  for  the  parameter,  it  can  be 
noticed  that  there  is  a  large  deviation  in  the  initial  part  of  the 
simulation.  This  deviation  is  caused  by  the  error  in  the  initial 
states.  Nevertheless,  the  on-line  identified  parameter  is  close 
to  the  off-line  determined  value.  Therefore,  the  performance 
of  DEKF  estimator  is  comparable  to  that  of  KF  estimator. 
In  particular,  the  RMSE  for  core  temperature  estimation  by 
DEKF  is  0.26,  the  same  RMSE  by  the  KF.  Even  though  there 
is  a  slight  error  between  the  measured  and  estimated  temper¬ 
ature,  the  closed  loop  estimators  show  a  good  performance  in 
predicting  core  temperature  overall.  As  discussed  in  Section 
III-C,  thermal  properties  can  vary  with  respect  to  operating 
temperature.  Therefore,  it  is  expected  that  better  performance 
can  be  achieved  by  using  temperature-dependent  parameters 
for  thermal  conductivity  and  specific  heat  coefficient. 

Figure  1 1  illustrates  the  performance  of  temperature  estima¬ 
tion  by  the  closed  estimator  in  stage  II  when  there  are  sudden 
changes  in  the  cooling  condition.  The  KF  can  accurately 
estimate  the  core  temperature  with  information  about  the 


change  in  parameter  value.  Since  the  DEKF  is  capable  of 
compensating  inaccuracy  in  the  parameter  of  the  system,  the 
DEKF  provides  reasonably  accurate  estimates  for  the  core 
temperature  by  comparing  the  core  temperature  predicted  by 
the  KF  estimator.  Even  though  the  RMSE  for  core  temperature 
estimation  by  DEKF  is  slightly  larger  than  the  RMSE  by  the 
KF,  the  error  is  still  reasonably  small  with  considerations  of 
the  sensor  accuracy. 

As  seen  from  Fig.  12,  the  KF  estimator  overestimates 
the  core  temperature  when  the  incorrect  parameter  value  is 
used  as  a  convection  coefficient.  In  other  words,  the  reliable 
estimation  of  core  temperature  is  only  possible  when  the 
accurate  parameter  is  available.  Thus,  it  can  be  concluded 
that  the  DEKF  estimator  outperforms  the  KF  estimator  due  to 
the  capability  of  parameter  identification.  The  RMSE  for  core 
temperature  estimation  in  stage  III  can  be  substantially  reduced 
from  1.18  to  0.31  by  the  DEKF.  It  is  worthy  to  note  that  the 
DEKF  can  be  augmented  with  other  existing  fault  detection 
methods  and  power  management  strategies  to  improve  the 
system  robustness  without  cost  increase.  For  instance,  in  order 
to  detect  partial  blockage  in  cooling  system,  typically,  a  mass 
flow  or  pressure  sensor  is  required.  The  DEKF  enables  the 
identification  of  convection  coefficient  by  using  sensors  which 
are  already  equipped  with  the  battery.  Thus,  by  monitoring  the 
difference  between  the  identified  and  off-line  predetermined 
values,  the  malfunction  of  cooling  system  can  be  detected. 

TABLE  m 

Performance  of  temperature  estimation:  RMSEs  for  core  and 

SURFACE 


Method 

DEKF 

KF 

Location 

Core 

Surface 

Core 

Surface 

Stage  I 

0.26 

0.07 

0.26 

0.07 

Stage  II 

0.39 

0.08 

0.29 

0.08 

Stage  III 

0.31 

0.11 

1.18 

0.15 
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This  performance  tells  that  a  significant  rise  in  temperature  can 
be  prevented  by  limiting  maximum  discharge/charge  current 
rate  computed  from  the  estimated  parameter  and  states.  Thus, 
the  operation  of  the  battery  system  can  be  safe  and  robust. 


VI.  CONCLUSION 

In  this  study,  a  method  to  estimate  the  temperature  distri¬ 
bution  in  cylindrical  batteries  under  unknown  cooling  con¬ 
dition  is  proposed.  First,  a  radially  distributed  1-D  thermal 
modeling  approach  for  a  cylindrical  battery  cell  is  considered 
and  polynomial  approximation  is  applied  to  obtain  a  reduced 
order  model  enabling  the  development  of  real-time  applica¬ 
tions.  Frequency  domain  analysis  shows  that  the  proposed 
model  provides  sufficiently  accurate  prediction  of  core  and 
surface  temperatures  with  a  reasonable  assumption  that  the 
temperature  of  cooling  media  does  not  change  rapidly.  The 
proposed  model  is  used  to  identify  thermal  properties  and 
convective  coefficient  for  a  2.3  Ah  26650  LFP  battery  cell 
using  a  set  of  measured  data.  The  numerical  analysis  on 
parameter  sensitivity  supports  the  use  of  constant  parameters 
for  thermal  conductivity  and  heat  capacity  and  the  impor¬ 
tance  of  identifying  the  convection  coefficient  on-line.  Then, 
the  Dual  Extended  Kalman  Filter  is  applied  to  estimate  the 
temperature  inside  the  battery  and  convection  coefficient  by 
the  cooling  fan.  The  proposed  method  requires  no  knowledge 
of  the  convective  cooling  conditions.  The  results  shows  that 
the  proposed  DEKF  estimator  can  provide  reasonably  accurate 
estimates  of  core  temperature  and  convection  coefficient  by 
using  surface  temperature  which  is  relatively  easy  to  measure 
in  practice.  In  addition,  a  faulty  operation  in  the  cooling  system 
can  be  detected  by  monitoring  the  difference  between  the 
identified  and  off-line  predetermined  values.  Since  forced  air 
is  used  as  a  cooling  media  to  reject  heat  from  the  cell  in  this 
paper,  the  range  of  the  convection  coefficient  of  which  we  are 
interested  is  less  than  100  W/m2K.  Therefore,  the  reader  is 
urged  to  investigate  whether  the  polynomial  approximation  is 
valid  for  their  applications. 

In  the  future,  the  proposed  method  can  be  used  to  develop 
various  battery  management  strategies,  e.g.  the  determination 
of  maximum  current  with  consideration  of  thermal  constraints 
or  optimal  fan  scheduling  for  energy  efficiency,  leading  to  the 
safe  and  efficient  operation  of  the  battery  system. 


Appendix  A 


The  matrices  A([k_1  ,Ck,  and  Ck  are  calculated  by 
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